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Abstract 

We derived an analytic formula for T violation by using the perturbation theory for 
small quantities, Am\^L/2E and 6a{x)L/2E, where 6a{x) represents symmetric and 
asymmetric matter fluctuations, i.e., deviations from the average density. We analyzed 
the effect of matter fluctuations to T violation, by assuming PREM profile of earth 
matter density. We found that matter fluctuations do not give any viable contribution for 
L <6000km, while the fluctuation effect becomes large due to resonances for L >7000km. 
For 7000km< L <8000km, matter fluctuations contribute destructively to the average 
density term and the net result is small, while for L >8000km, the contribution from 
matter fluctuations becomes large but contributes constructively. 
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1 Introduction 



In our previous paper [1], we analyzed the matter fluctuation effect to T violation, P{i'a 
Vjj) — P{v0 — > Va) {oi. 7^ /3) at a neutrino factory[2] by using the perturbation method 
developed by Koike and Sato [3] and Ota and Sato[4]. The perturbation is made with 
respect to small quantities, ^m\iL/2E and 6a{x)L/2E, where 6a{x) represents matter 
fluctuations, i.e., deviations from the average density . We examined T violation up to 
the 2nd order and found that it arises from the average density, the 1st order term which 
is proportional to A'mliL/2E and term from matter fluctuations, the 2nd order term 
proportional to {Am2iL/2E)(da{x)L/2E). The zeroth order term and terms proportional 
to {Aml^L/2Ef and {6a{x)L/2E)'^ (n = 1,2,3) do not contribute to T violation. The 
1st order term and the 2nd order term from symmetric matter fluctuations which we 
denote Sag contribute to sin 5 term and the 2nd order term from asymmetric matter 
fluctuations, 5aa does to the fake cos 6 term, where S is the CP violation phase in MNS 
neutrino mixing matrix [5]. 

By using the preliminary reference earth model (PREM) [6] for symmetric matter fluc- 
tuations and assuming that asymmetric matter fluctuations are much less than symmetric 
matter fluctuations given by PREM, we computed T violation and found that the 2nd 
order term from symmetric and asymmetric matter fluctuations gives only negligible con- 
tributions to T violation, and thus the constant (average) matter approximation is valid 
for L = 3000km. On the other hand, for L = 7332km, we found that the contribution 
from symmetric matter fluctuations becomes as large as the 1st order term and moreover 
they contribute destructively so that T violation becomes very small. This means that 
the constant (average) matter approximation is not valid for L >7000km and also the 
validity of our 2nd order formula should be examined. 

In this paper, we discuss the following three questions: (1) Is the 2nd order formula 
valid? (2) What is the length where the constant (average) matter approximation fails for 
T violation? That is, at what length, the matter fluctuation effect becomes important. 
(3) What is the size of T violation for L >7000km? 
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To answer these questions, we computed the next order (3rd order) contribution to 
T violation, i.e., the term proportional to (Aml-^^L /2E) (6a(x) sL /2Ey . Our result is as 
follows: The contribution from matter fluctuation can be safely neglected for L <6000km. 
When we discuss T violation with length larger than 6000km, the matter fluctuation effect 
should be taken into account. The 3rd order contribution is negligible in comparison 
with the 1st and the 2nd order term for all distances. Therefore, T violation can be 
rehably estimated for all distances and it becomes very small for 7000km< L < 8000km. 
For L >8000km, the 1st and the 2nd terms contribute constructively and T violation 
becomes large. 

This paper is organized as follows: In Sec. 2, the analytic formula for T violation is 
given. The numerical analysis for T violation by using PREM profile is given in Sec.3 
and the mechanism how the cancellation occurs for 7000km< L < 8000km is explained 
analytically. The summary is given in Sec. 4. The derivation of the analytic formula is 
given in Appendix. 



2 T violation formula 



For completeness, we give the 1st and the 2nd order contributions to T violation and the 
definition of parameters in the formula. We also give the 3rd order contribution from 
symmetric fluctuations, 
(a) Notation 

We begin with defining the neutrino mixing matrix as 



U = e^^^^Miag(l, 1, e^'')e^^^^^e^^-^^ 



-SxCy CxSyS^G 



CxCy SxSyS^C 



i5 ^ ^ Ab 



(1) 



where \ {j — 2,5,7) are Gell-Mann matrices and Ca = cos 9 a and Sa — sin^^. Since 
the Majorana CP-violation phases are irrelevant to the neutrino oscillations (flavor os- 
cillations) [7] , we neglected them. The relation between the flavor eigenstates, I^q) {a — 



e, II, t), and the mass eigenstates, ji/j) {i — 1, 2, 3), is given by 

Wa) = Uai\iyi) . 

The evolution of the flavor eigenstates in matter with energy E is given by 
where Hamiltonian H[x)i3a is given by 



(2) 
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Here Am?- = mf —ni^ with mj being the mass of ji/j), is the Fermi couphng constant 
and 

a(.) . 2v^G.„.(.)E ^ 7.56 x 10"= (^) (^) eV^ , (5) 

where ne{x), Yg and are the electron number density, the electron fraction and the 
matter density, respectively For the electron fraction, we use Y^ — 0.5. 

We separate the matter density fluctuation from its average a, 



5a{x) = a{x) — a , 



(6) 



and consider the deviation 5a{x) as a perturbative term. That is, we solve the evolution 
equation by treating 6a{x)L/2E and /^m\iL/2E as perturbative terms, because they are 
small for most of the cases of planned neutrino factories. 
T violation is deflned by 



AP^ 



(7) 



which is evaluated by using the method developed by Koike and Sato [3], Ota and Sato [4]. 
Ota and Sato showed that the 1st order approximation for da{x)L/2E is good enough 
to reproduce the transition probability. However, T violation is as small as a few % of 
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the transition probability so that this approximation is not vahd. In fact, we showed the 
2nd order term becomes as important as the 1st order term for L > 6000km. We show 
T violation, AP^^^^^ up to the 3rd order perturbation for symmetric matter fluctuations 
and the 2nd order for asymmetric matter fluctuations. The analytical formula is given 
by expanding symmetric and asymmetric matter fluctuations in terms of Fourier series 
as 

5a{x)s^ E «2ne-'^^"^ 5a{x)a^ E «2n+le-'^^"+^^ (8) 

n=±l,±2,- n=0,±l,- 

where 



riTT 



Qn = 



(9) 



(b) The analytic formula for T violation 

In order to define T violation, we define the following quantities: 



tan 29z 
A± 

k 



C2z{Amli - Amlisl) - a ' 
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2E 
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2E 
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2E 



(10) 



The sum of the 1st and the 2nd order terms due to symmetric matter fluctuations 
contributes to sin 5 term and is given by 



A pT{l+2s) 
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0'2n 
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(11) 
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The asymmetric matter fluctuation contributes to the cos 5 term and is given by 



A pr(2a) 



X 
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E 



+ 



kiL k2L kL 
sm — — sm — — sm — 
2 2 2 



k 



n=0,l,- 



cot — - 
2E k'^ — q2n+i 2 k\ — 



cot 



k^L 
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U2 _ „2 



cot 



koL 



^2 ?2n+l 2 J 

The 3rd order contribution from symmetric fluctuations is given by 



(12) 
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(13) 



It is amusing to see that the 1st and the 2nd order terms due to symmetric matter 
fluctuations have the same coefficient and the same oscillation term as we can see from 
Eq.(ll). If one of resonance conditions, kf — A;| — g| and k"^ — g| is realized at 
some distance, the matter fluctuation term dominates over the 1st order term, although 
this singular behavior is cancelled by the oscillation term. The asymmetric matter con- 
tribution gives the similar contribution except for the angle 5. The 3rd order term is 
rather complicated and seems to have double poles when n = m, which is false due to 
the cancellation between other terms, and there are no singularities. 

As we expected, symmetric and asymmetric matter fluctuations contribute to the 
sin S and cos 5 parts of T violation. 

The important fact for T violation is that the relation 



AP^ = AP^ = AP^ 



(14) 



holds even in our 3rd order formula. That is, we can not gain any further information by 
examining other channels. This fact for the constant matter was first found by Krastev 
and Petcov[8] and the validity for the 2nd order formula is proved in Ref.l. 

3 Numerical analysis 

By using the analytic formula, APj^[,^^+^*+^"+^*\ we investigate the L and E dependences 
of T violation. The energy of neutrino, E and the distance, L are taken from IGeV to 
30GeV and from 1000km to 12000km. Neutrino oscillation parameters are taken as 

Am^i = 3 X 10-2 (cV^) , Aml^ = 5 x 10"^ (eV^) , 
sin2^^ = sin2^2/ = l' sin^^-0.1, 5 = 7r/4 , (15) 

as a typical case. 

(a) General features of T violation 

Here we examine the effect of symmetric fluctuations to T violation. For the average 
density and symmetric fluctuations, we assume PREM and decompose it into the average 
density a and the Fourier coefficients a2n {p, = 1,2,3,4), which are functions of L. It is 
a general belief that the average density is determined with better accuracy than matter 
fluctuations. 

We computed T violation, APj'^^ from L =1000km to 12000km. We found that the 
2nd and the 3rd order terms are not important for L < 6000km and thus the constant 
(average) matter approximation is valid. We remind that the 1st order term contains the 
effect of the constant matter, while the 2nd and the 3rd order terms are due to matter 
fluctuations as we can see in Eqs.(ll), (12) and (13). 

Now we concentrate on the distance L > 6000km. In Figs. 1-6, we show the E depen- 
dence of T violation for L =6000km, 7000km, 7700km, 8000km, 10000km and 11000km, 
where the 1st, the 2nd and the 3rd order terms are shown by the dashed line, the dotted 
hue and the dash-dotted hue, respectively. The sohd hue represents the sum of them. 
The vacuum (no matter) case is shown by the dash-twodotted line for comparison. 
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When L > 7000km, the 2nd order term becomes the same order as the 1st order term 
and is no more neglected. Moreover, as we see in Figs. 2 and 3, for 7000km< L <8000km, 
the 1st and the 2nd order terms contribute destructively and T violation becomes small. 
When 8000km< L <9000km, the 1st and the 2nd order terms contribute constructively 
as we see in Fig.4. For L > 10000km, the 2nd order term becomes much larger than the 
1st order term. 

For all distances, the 3rd order term (the dash-dotted line) gives only a negligible 
contribution and it is hard to see from these figures. Thus, the 3rd order term can safely 
be neglected so that we need to consider the 1st and the 2nd order terms only. 

When we see Figs. 1-6, for E > 5GeV, we observe that there are two peaks. One is at 
E =5.5GeV and the other is at E =10GeV for L =6000km. The energies of these peaks 
increase as L increases. T violation for energies larger than that of the lower energy peak 
lies between these two peak values. 

We examined the L dependence of these two peak values and the result is shown in 
Fig. 7a. Faint diamonds and dark diamonds correspond to the peak values of the lower 
and the higher energy peaks, respectively. Diamonds, stars, boxes, triangles, circle show 
T violation for = sin^ig = 0.1, 0.08, 0.06, 0.04 and 0.02. Fig.7b shows the case for 
E =30GeV. 

Due to the cancellation between the 1st and the 2nd order terms, T violation becomes 
zero at around 7300km. For E =30GeV, there are two zeros at around L =7300km and 
9600km. 

Another aspect we can observe from these figures is that the dependences of T 
violation is roughly linear. 

(b) The relative sign between the 1st and the 2nd order terms 

In order to see the changes of the relative sign between the 1st and the 2nd order 
terms, we write these terms omitting the overall factor and the oscillation term as 



The 1st order term is expressed by 1, while the 2nd order term is expressed by the 




(16) 
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sum of Fourier coefficients, a2n- Tlie singularities at — A;f = and k'^ — 
correspond to resonances. In general, the 2nd order term is smaller than the 1st order 
term if ki, ^2 and k arc away from resonance points. Since q2n = 2mr/L, we compared 
{kiL/2y, {k2L/2f and {kL/2f with {m:f for E >5GeV in Fig.8. For L is as small 
as 3000km, the resonance condition is not satisfied so that the 2nd term gives only a 
negligible contribution. 

We consider the change of the relative sign between the 1st and the 2nd order terms 
as L increases from 7000km to 8000km at the lower energy peak position. Since 02 < 0, 
/ci < 0, A;2 > and A; > 0, the sign of the ki term is negative for L= 7000km and becomes 
positive at L =8000km. The resonance occurs at around L =7700km. The same is true 
for k2 term. The k term never reaches to the resonance point. That is, when L is as 
small as 6000km, all terms can be neglected because they never approach to the resonance 
point. As a result, T violation is positive. When L exceeds 7000km, contributions from 
the ki and /c2 terms become important and their signs are negative. As L approaches 
to 7700km, their contribution cancelles the 1st order term and T violation vanishes at 
around 7550km. Between 7550km< L <7700km, both ki and k2 terms become still 
negative and dominates over the 1st term so that T violation becomes negative. At 
the resonance point, this singularity is cancelled by the oscillation term and T violation 
obtains non-zero negative value. After passing the resonance point, both ki and k2 terms 
become positive and contribute additively to the 1st order term. Since the oscillation 
term becomes negative and thus T violation remains to be negative. 

The L dependence of T violation at the larger energy can be understood similarly. 
For E =30GeV, the zero of T violation at around 7300km can be understood similarly. 
The zero for 9600km is due to the resonance for k. 

(c)The effect from the uncertainty for the average matter density 

So far, we used PREM to derive the average density and the matter fluctuations. For 
the sake of argument, we consider 5% uncertainty for the average matter density although 
we do expect that the uncertainty is much less. We examined how T violation changes 
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when the average density is changed by 5%. If the average matter density changes, it 
affects to the angle Sz and ki, /c2 and k as we can see in Eqs.(lO). As a result, the 
distance L where the resonance occurs changes. In Fig.9, we show the L dependence of 
T violation. Diamonds, boxes and stars correspond to the average value from PREM, 
5% smaller value and 5% larger value. Faint and dark ones correspond the lower energy 
peak and the higher energy peak. It is interesting to see that the zero point shift to the 
longer (shorter) L by about 200km as the average density becomes smaller (larger). 

(d) The effect from the uncertainty of matter fluctuations 

As we discussed in the subsection (b) , the main contribution from matter fluctuations 
is from the term containing the Fourier coefficient 02 for symmetric fluctuations, and 
similarly ai for asymmetric fluctuations. Since 02 is determined from the most dense 
matter part (middle part) along the neutrino path. On the other hand, higher modes a2n 
(n = 2, 3, • • •) are determined mainly by the crust of earth. Therefore, for L > 6000km, 02 
is considered to be rather unambiguously determined reflecting the deep inside structure 
of mantle. 

For asymmetric fluctuations, Oi reflects the global asymmetric feature of matter pro- 
file. For distances at a neutrino factory, neutrinos pass mainly through the mantle and 
we do not expect much uncertainty ai. For the shorter length, neutrinos pass through the 
crust and sometimes the sea. In this situation, the matter profile with large asymmetric 
matter fluctuations may need to be considered. [9] 

For the sake of argument, we assume that there is about 10% uncertainty for 02, i.e., 
Q2n = (1 =t ^■^){0'2n)pREM- Sincc APj^^^^ depends linearly on 02, the 10% uncertainty 
for 02 gives the same uncertainty for APjj^^\ For distances where the 2nd order term 
is neglected, the uncertainty from symmetric fluctuations is negligible (L <6000km). In 
distances where the 2nd order term dominates, then 10% uncertainty appears (L > 
10000km). The uncertainty for L > 8000km is smaller than 10%. For 7000km< 
L <8000km, the uncertainty is larger than 10%. In Fig. 10, we show how T violation 
becomes uncertain if 02 has 10% uncertainty. 
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For asymmetric fluctuations, we assumed that a2n-i = 0.1{a2n)pREM in addition 
to symmetric matter fluctuations determined from PREM. The effect from asymmetric 
fluctuations is very similar to the case of symmetric fluctuations and it is shown in Fig. 10. 

4 Summary 

In this paper, we derived the analytic formula for T violation up to the 3rd order term 
for small quantities, /\m\^L/2E and 6a{x)L/2E. By using this formula, we examined 
the E and L dependence of T violation. 

We showed that the effect from both symmetric and asymmetric matter fluctuations 
are negligible for L < 6000km and the constant (average) matter approximation (the 1st 
order term) is valid. Therefore, T violation contains the uncertainty from the average 
matter aside from mixing angles. Since the average is considered to be determined with 
much less uncertainty than matter fluctuations, T violation is determined with a good 
accuracy for L < 6000km. 

For L > 6000km, the situation changes. Matter fluctuations (the 2nd order term) 
give a sizable effect to T violation. Moreover, the 1st and the 2nd order terms contribute 
destructively for 7000km< L <7700km. As a result, the T violation becomes very small. 

In Fig. 7, we showed the L dependence of T violation at the lower energy peak 
(the peak at E =5.5GeV for L =6000km) and the higher energy peak (the peak for 
E =10GeV for L =6000km). T violation at the lower energy peak has the largest value 
for 5000km< L <6000km and L ~10000km and T violation at the higher energy peak 
docs for 3000km< L <4000km. T violation becomes zero at around L =7300km. This is 
due to the resonance effect, which we explained why this happens by using the analytic 
formula. For E — 30GeV, T violation behaves similarly and has the largest value at 
5000km. 

We also examined how the L dependence of T violation varies as the average matter 
is changed by ±5%, although we believe that the average matter density is determined 
with much less uncertainty. We found that the distance which gives zero of T violation 
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shifts about ^200km. 
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Appendix A : The brief summary of our previous re- 
sult 

In Rcf.l, wc derived T violation fromula up to the 2nd order with respect to small 
quantities of Am|iL/4£^ and SaL/4E, where Sa represents matter fluctuations, where 
both symmetric Sus and asymmetric Sua fluctuations are taken into account. Here, we 
give the brief summary of the previous result and the 3rd order calculation needed to 
compute the next order symmetric matter fluctuation effect. 

We express the S-matrix as 

S = Soo+ E Sj-f, (A.l) 

n,m=0,l,2,--- 

where 6*00 is the 0th order term from Hqq, SqI'"^^ represents the (n+m)th order term from 
H^^ X H^, i.e., the term of order (Am^iL/4£;)"((5aL/4£;)"^. 

The contributions to T violation is summarized as follows: (1) The 1st order term 
is from SooSq\'^* and is proportional to sin 5 and contains the contribution from the 
constant (average) matter. (2) The 2nd order term is from S'oo'S'oi'i''* and 5'oiiS'o*]''^''*. 
The symmetric and asymmetric matter fluctuations contribute to sin 5 and cos 5 terms, 
respectively. (3) The 2nd order term, and the 3rd order term, do not 

contribute to T violation. 

Firstly, we present the result given in Ref.l: 

C250- SyS2ze'~''^ (})- CyS2ze''''^(f)- \ 

SyS2ze'^(f)- 0+ + Slc2z(f>- - C^(02- " 01-) S2y (c|02- " s|0i_) 

\ CyS2ze''^(p- S2y (ci02- - Si01-) <P+ + CyC2i0_ - ^(^2- " 01-) / 

(A.2) 



2 



where 0± and 0i± {i — 1,2) are deflned as follows. 

0i± = e-'"""^ ± e-'"-^ , 02± = e~^"+-^ ± e"'""-^ . (A.3) 
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5'oi'i"^(n > 1) is parametrized by 
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(A.ll) 



By using the above fromula, we obtained the 2nd order formula for T violation. 

Finally, we give the result for 5'q?'i^ which is needed to estimate the 3rd order correc- 

:.(o,2) 



tion. 5'oi'i is calculated from 



q{0,2) p(0,2) ^ 



(3(0,2) _ p(0,2) 



S2zC2z 

{2Ef 



E 



\in+n=even 



{k + qn){k + qn + qm). 



S25C25 



E 



{2Ey I +t^„dd + <ln){k + qn + qm) 



14 



£:(o,2) _ jr(0,2) 



'25 



m 



((-!)" + l)((-ir + l)a„a, 
{k + qn){k + qm) 



(A.12) 



Appendix B : The 3rd order calculation 

We consider the 3rd order correction form HqiH^ where we take symmetric matter 
fluctuations for Hi. The contribution arises from 5'oo5'oi'i'**, Sq\'^i SqI'^* and 5'qi'}'*S'o°'i^*. 
The latter two can be calculated the result given in Appendix A. Here, we give the result 
for Sqi i which is given from 

92m ^2 ~ 92m J ki + 92n+2m / 



k, (5 (2E)2 A: + g^n 



e 



»(l,2s) _ B/(l,2s) _ _Arn2iS2x^ 



X 



2 2 



/^i A;2 J [rT'o (2£^)^ V^i + 92m A;2 - g2m J k2- q2n 



SzSz-z 2t I 0'2nO'-2n 1 



+2m 



-c: 



1^ E 



(2^)2 k + q2n 



(B.l) 

By using the above formula, we can derive the 3rd order contribution to T violation. 
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Figure 1: The energy dependence of the T violation, AP^^^ with L — 6000km for 
(a) 5GeV< E <30GeV and (b) lGeV< E <5GeV. The dashed hne, the dotted hne 
and the dash-dotted hne show the 1st, the 2nd and the 3rd order terms. The sohd 

hne represents the sum of them. The vacuum (no matter) case is shown by the dash- 
twodotted hne for the comparison. In this plot, we use sin2^j; — sm20y = 1, sin^^ — 0.1, 
Amii = 5 • lO-^eV^ Ami^ = 3 • lO'^eV^ and 5 = 7r/4. 
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Figure 2: The energy dependence of the T violation, AP^^,^ with L — 7000km for (a) 
5GeV< E <30GeV and (b) lGeV< E <5GeV. The species of the hnes and the oscillation 
parameters are the same as those in Fig.l. 
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Figure 3: The energy dependence of the T violation, AP^^,^ with L — 7700km for (a) 
5GeV< E <30GeV and (b) lGeV< E <5GeV. The species of the hnes and the oscillation 
parameters are the same as those in Fig.l. 
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Figure 4: The energy dependence of the T violation, AP^^,^ with L — 8000km for (a) 
5GeV< E <30GeV and (b) lGeV< E <5GeV. The species of the hnes and the oscillation 
parameters are the same as those in Fig.l. 
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Figure 5: The energy dependence of the T violation, AP^^,^ with L — 10000km for (a) 
5GeV< E <30GeV and (b) lGeV< E <5GeV. The species of the hnes and the oscillation 
parameters for (a) 5GeV< E <30GeV and (b) lGeV< E <5GeVare the same as those 
in Fig.l. 
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Figure 6: The energy dependence of the T violation, AP^^,^ with L = 11000km for (a) 
5GeV< E <30GeV and (b) lGeV< E <5GeV. The species of the hnes and the oscillation 
parameters are the same as those in Fig.l. 
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Figure 7: The L dependence of T violation, A^ ,^^. In Fig. (a) T violation at the lower 
energy peak {E = 5 ^ 8 GeV) which is shown by faint points and the higher energy peak 
[E = 10 ^ 13 GeV) which is shown by dark points. Diamonds, stars, boxes, triangles 
and circles show T violation for — sin ^13 — 0.1, 0.08, 0.06, 0.04 and 0.02 respectively. 
In Fig. (b) we show the values of T violation AP^^^^ at = 30 GeV. The oscillation 
parameters except for Sz are the same as those in Fig.l. 
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Figure 8: The energy dependence of (A;iL/2)2(Fig.(a)), (A;2V2)^(Fig.(b)) and 
(fcL/2)^(Fig.(c)). Dotted, dashed, dash-dotted and dash-twodotted hnes correspond to 
L = 3000,7000,7700,11000 km, respectively. The sohd line shows tt^. The oscillation 
parameters are the same as those in Fig.l. 
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Figure 9: The average density (a) dependence of T violation, AP^^,^. Diamonds, boxes 
and stars correspond to the average density value from PREM, 5% smaller value and 5% 
larger value. Faint and dark points correspond to the lower energy peak (E — 5 8 
GeV) and the higher energy peak [E = 10 ^ 13 GeV). The oscillation parameters are 
the same as those in Fig.l. 
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Figure 10: The variation of T violation, AP^^^ by the deformation of the shape of 
matter profile for L = 7700km. Dotted lines show changes of T violation when Fourier 
coefficients for symmetric matter fluctuations are varied by 10% from PREM. Dash- 
dotted lines are those when asymmetric matter fluctuations are associated to be about 
10% of PREM symmetric matter fluctuations. The oscillation parameters are the same 
as those in Fig.l. 
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